

IMAGE RPI Reawakens Plasmaspheric Refilling Research 
D. L. Gallagher and Z. B. Smith 

The plasmasphere is a toroidal region of cold plasma surrounding the Earth that results 
from ionospheric outflow and accumulation. The physics of refilling and the dynamics of 
this region have been studied for nearly 50-years. During that time many models have 
been proposed, but little has been done to test these models due to a lack of observational 
information. With the launch of the IMAGE Mission in March 2000 the Radio Plasma 
Imager has provided true field aligned density measurements that uniquely enable the 
testing of these models and a final determination of the physical processes important for 
the plasmasphere ’s recovery from storm-time conditions. 
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GLOBAL DYNAMICS OF THE EARTH’S PLASMASPHERE 
PhD Dissertation, Maria Spasojevic, 2003. 
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An Abridged History of Plasmaspheric Refilling 
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The Birth of Refilling 

(first there came emptying) 


The need to understand 
plasmaspheric refilling can be traced 
to the first observations of 
plasmaspheric erosion. Using 
charged particle traps on the Soviet 
lunar vehicles Luna 1 and Luna 2, 
Gringauz [1959] made the first in 
situ measurements of thermal ions in 
the “outermost ionosphere” and 
using ground whistler measurements 
Carpenter [1963] discovered an 
outer boundary to this thermal 
plasma region found to be sensitive 
to geomagnetic activity. 


N 

(•l/cc) 



(b) 

Overall loss or redistribution? 


How Will Plasma 
Become Distributed 
Along Field Lines? 


Taylor et al., [1965] concludes their Ogo-A 
ion measurements are consistent with the 
models of field aligned density by Johnson 
[1960], Bates and Patterson [1961], and 
Angerami and Thomas [1964] at 2000 km, 
but “exhibit significant departure from 
these models” at high altitudes up to 30,000 
km. 

Equatorial densities don’t fit the 
Diffusive Equilibrium model at 
high altitude. 
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Ionosphere - Plasmasphere Exchange 


Parks [1970] found that 
daytime ionospheric outflow 
was more than enough to 
allow down flow from the 
plasmasphere at night to 
maintain ionospheric 
densities. Flux tubes 3.5<L<5 
take 8-days or more to refill, 
but probably never get there 
due to resumed activity. 
Plasmaspheric recovery 
comes from ionospheric 
outflow. 
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The ionosphere is the source for (slow!) plasmaspheric refilling. 


Hydrodynamic Top-Down Refilling 



The single-stream hydrodynamic model of Banks [1971] shows the evolution of 
supersonic streams from the ionosphere during the early stages of refilling to 
subsonic flows and eventually to a diffusive equilibrium. In the initial phase, 
supersonic flows from conjugate ionospheres interpenetrate at the equator of the 
depleted flux tube causing the formation of shocks. The second stage subsonic 
refilling occurs once the shocks reach the topside of the ionosphere in conjugate 
hemispheres. 

Interacting interhemispheric shocks: Fast two-stage refilling 
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Two-Stream Top-Down 



Two-stream models show or do not show shock formation, dependent 
on the physics of the inter-penetrating interhemispheric streams. 


Interacting interhemispheric shocks: Slower refilling 


Two Stream Anisotropic 
Top-Down Refilling 


Singh and Torr in 1990 proposed a two stream 
anisotropic (TAN) model with equatorial ion 
heating by wave-particle interaction. TAN 
models show a two-stage refilling. The first 
stage is characterized by development of 
supersonic flows in depleted flux tubes, 
interhemispheric plasma exchange, shock 
formation in each plasma stream and rapid 
propagation of the shock first upward to the 
equator and then downward to the ionosphere. 
This first stage lasts a short time (~5 hrs) with 
minimal refilling taking place. When the shocks 
reach the ionosphere, a quasi-steady distribution 
results in the flux tube. This allows the second 
stage of refilling from the ionosphere to occur. 
These flows are generally subsonic and isotropic 
in nature. 
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Non-interacting interhemispheric shocks with slow refilling. 


Density (cm 3 ) 


Semi-Kinetic Bottom-Up Refilling 



In 1992, Wilson used a semi-kinetic approach that includes 
Coulomb collisions and wave-particle interactions to model 
early-time plasmaspheric refilling. 


Field aligned profile fills without shock structure. 



Density Variation Along the Magnetic Field 
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We could only approximate the field-aligned density profile. 
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True Field- Aligned Densities 



Physics of plasmaspheric refilling is enabled by RPI 

Obvious, but not anticipated! 
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First true measurements of field aligned density profiles 


Plasmaspheric Refilling Study 

Renewed 


Reinisch et al., [2003] apply, for the first 
time, RPI field aligned derived densities 
to the pattern of refilling. 
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Solar, Season, & Activity 
Dependent Refilling Rates 


Inspired by IMAGE RPI and EUV like 
many researchers, Su et al., [2001] use 
geosynchronous observations to find the 
refilling rate to be non-linear in time and 
dependent on geomagnetic activity, solar 
conditions, and seasons. 



K p 


Electron Dens ty (cm 


Activity Dependent 
Field-Aligned Density Gradients 


05 Novemoer 2004 07 Ncvemoer 2004 22 March 2001 



Tu el al., 2006 demonstrate the initial steps for the formation of an empirical model 
for plasmaspheric and magnetospheric trough field line densities as a function of 
geomagnetic conditions. 


New step toward field-aligned density empirical modeling. 



Field-Aligned Densities in 2001 
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Smith has derived 1556 field-aligned densities profiles from 
soundings in 2001 as a resource to study plasmaspheric refilling. 
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2-D Drift Velocities from the IMAGE 


EUV Plasmaspheric Imager 


D. Gallagher and M. Adrian 



Will focus discussion of the 
technique and initial results 
on one plasmaspheric 
erosion time period 


BEFORE AFTER 

05:06 08:00 



10 July 2000 (a) 





Regional Plasmaspheric Drift 

Objective: 

Use all available intensity structure in a series of EUV images to quantify movement of 
image elements within the field of view to fulfill promise of obtaining 2D E-fields. 

Preparation: 

Perform the analysis using image observations projected into the magnetic equatorial plane 
such that image products can be directly compared. That is done here using the pseudo- 
density calculation. 




To Compute 2D Piasmaspheric Drift 


Techniques: Optical Flow and Cross-Correlation 


Optical Flow: 


Optical flow constraint equation 
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Tried by Christophe Bernard (Ecole Polytechnique in France) 
without useful results so far. IDL code available at 
http://csrsrvl .fVnu.ucl.ac.be/csr web/data/ 


Correlation 
Array \ 


Cross-Correlation: 


With this technique image elements 
are directly compared for 
correspondence. 

Selected elements from one image or 
kernels are swept across portions of a 
second image looking for the best 
match. 

The shifted position of the best match 
defines the motion of that image 
element. 
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Here’s the Challenge 


• Plasmapause 

• Noise 

• Shadow 

• Ionosphere 

• Other features? 




Subarrays being compared 




Image for region selection 
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Noise Dominates - 

Without treatment, noise —4. 

dominates derived flow velocities 
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A Few Noise Removal Techniques Tried 
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Solar Wind E-Field 10 July 2000 





Initial Results for 
July 10, 2000 
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Analysis performed 
between image products 
separated by 20 minutes 

Noise removed from linear 
pseudo-density. 
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Correlation performed on 
linear pseudo-density. 
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Solar Wind E-Field 10 July 2000 
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Solar Wind E-Field 10 July 2000 
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Solar Wind E-Field 10 July 2000 
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Solar Wind E-Field 10 July 2000 
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Solar Wind E-Field 10 July 2000 
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Summary 



Cross-Correlation analysis appears 
capable of deriving convective flows 
broadly throughout the inner 
magnetosphere 

As suspected, even “simple” 
plasmaspheric erosion is not so simple 

The technique requires refinement prior 
to zealous use 

Zealous use of a proven technique is 
warranted 


